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For a family of logistic-like maps, we investigate the rate of convergence to the critical attractor 
when an ensemble of initial conditions is uniformly spread over the entire phase space. We found 
that the phase space volume occupied by the ensemble W(t) depicts a power-law decay with 
log-periodic oscillations reflecting the multifractal character of the critical attractor. We explore 
the parametric dependence of the power-law exponent and the amplitude of the log-periodic 
oscillations with the attractor's fractal dimension governed by the inflexion of the map near 
its extremal point. Further, we investigate the temporal evolution of W(t) for the circle map 
whose critical attractor is dense. In this case, we found W(t) to exhibit a rich pattern with a slow 
logarithmic decay of the lower bounds. These results are discussed in the context of non-extensive 
I , Tsallis entropies. 
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. I. INTRODUCTION 

CO 

Nonlinear low-dimensional dissipative maps can describe a great variety of systems with few degrees of freedom. 
The underlying nonlinearity can induce the system to exhibit a complex behavior with quite structured paths 
in the phase space. The sensitivity to initial conditions is a relevant aspect associated to the structure of the 
dynamical attractor. In general, the sensitivity is measured as the effect of any uncertainty on the system's 
variables. For systems exhibiting periodic or chaotic orbits, the effect of any uncertainty on initial conditions 
depicts an exponential temporal evolution with = lim^^o)— >o Ax(i)/A:r(0) ~ e At , where A is the Lyapunov 
exponent, Ax(0) and Ax(t) are the uncertainties at times and t. When the Lyapunov exponent A < 0, 
characterizes the rate of contraction towards periodic orbits. On the other hand, for A > 0, it characterizes the 
rate of divergence of chaotic orbits. At bifurcation and critical points (i.e., onset to chaos) the Lyapunov exponent 
A vanishes. Recently, it was shown that this feature is related to a power-law sensitivity to initial conditions on 



the form |-§] 



m = ii + (i-qM 1/(1 - q) , (1) 



with X q defining a characteristic time scale after which the power-law behavior sets up. 

A quantitative way to measure the sensitivity to initial conditions is to follow, from a particular partition of the 
phase space, the temporal evolution of the number of cells W(t) occupied by an ensemble of identical copies of the 
system. For periodic and chaotic orbits W(t) — W(Q)e xt . In the particular case of equiprobability, the well-known 
Pesin equality reads i4T = AifA>0|| with K being the Kolmogorov-Sinai entropy Q defined as the variation 
per unit time of the standard Boltzmann-Gibbs entropy S — — ^2p% lnp;. This equality provides a link between 
the sensitivity to initial conditions and the dynamic evolution of the relevant entropy. 
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At bifurcation and critical points and for an ensemble of initial conditions concentrated in a single cell, i.e. 
W(0) = 1, it has been shown that 



W(t) = [l + (l-q)K q t} 1 /^ , (2) 

with K q being the generalized Kolmogorov-Sinai entropy defined as the rate of variation of the non-extensive 
Tsallis entropy S q = (1 - - 1) §• The Pesin equality can be generalized as K q = \ q if X q > M. 

Tsallis entropies have been successfully applied to recent studies of a series of non-extensive systems and provided 
a theoretical background to the understanding of some of their unusual physical properties |t],|8| . 

The expansion towards the critical attractor of an ensemble of initial conditions concentrated around the inflexion 
point of the map can be characterized by a proper S q evolving at a constant rate. Scaling arguments have shown 
that the appropriate cntropic index q is related to the multifractal structure of the critical dynamical attractor by 



(3) 



are the extremal singularity strengths of the multifractal spectrum of the critical attractor 
j. The above scaling relation has been shown to hold for the families of generalized logistic and circle maps 

HID- 



However, the temporal evolution of critical dynamical systems can be strongly dependent on the particular initial 
ensemble. Although some scaling laws can be found for an ensemble of initial conditions concentrated around the 
map inflexion point, these are usually not universal with respect to a general ensemble. In this work, we are going 
to numerically investigate the critical temporal evolution of the volume of the phase space occupied by an ensemble 
of initial conditions spread over the entire phase space. This ensemble is expected to contract towards the critical 
attractor. Using a family of one-dimensional generalized logistic maps having df < 1, we will perform a detailed 
study of the parametric dependence of W(t) on the fractal dimension of the critical attractor. Due to the discrete 
scale invariance of the critical attractor, the convergence displays log-periodic oscillations fDifl . We are also going 
to explore the dependence of the amplitude of these oscillations with respect to the attractor 's fractal dimension. 
Further, the behavior of W(t) will be investigated for the one-dimensional critical circle map having df = 1. For 
this map, the temporal evolution is expected to display distinct trends since the critical attractor is dense. 



II. THE CONVERGENCE TO THE CRITICAL ATTRACTOR OF GENERALIZED LOGISTIC MAPS 

Logistic-like maps are the simplest one-dimensional nonlinear dynamical systems which allow a close investigation 
of a series of critical exponents related to the onset of chaotic orbits. This family reads 

x t+ i = 1 - a\x t \ z ; (z>l ; 0<a<2 ;t = 0,1,2,... ;aj t g[-l,l]) . (4) 

Here z is the inflexion of the map in the neighborhood of the extremal point x = 0. These maps are well known to 
have topological properties not dependent of z. However, the metrical properties, such as Feigenbaum exponents 
p4] , p^| and the multifractal spectrum of the critical attractor do depend on z. In particular the fractal dimension 
of the critical attractor df(z) < 1 |l(| and therefore it does not fill a finite fraction of the phase space. For a 
set of initial conditions spread in the vicinity of the inflexion point, it was found that the volume in phase space 
occupied by the ensemble grows following a rich pattern with the upper bounds W max (t) governed by a power-law 
W max (t) oc t 1 /' 1 " 9 ', where q is the entropic index characterizing the relevant Tsallis entropy that grows at a 
constant rate. It has been shown that the dynamic exponent 1/(1 — q) is directly related to geometric scaling 
exponents related to the extremal sets of the dynamic attractor ||. 

Due to the presence of long-range spatial and temporal correlations at criticality, one expects the critical exponent 
governing the temporal evolution to be sensitive to the particular initial ensemble. Indeed, the multifractal 
spectrum characterizing the critical dynamical attractor indicates that an infinite set of exponents are needed to 
fully characterize the scaling behavior. In particular, an ensemble consisting of a set of identical systems whose 
initial conditions is spread over the entire phase-space is a common one when studying non-linear as well as 
thermodynamical systems. 

Here, we will follow the dynamic evolution, in phase space, of an ensemble of initial conditions uniformly 
distributed over the phase-space and explore its relation with the generalized fractal dimensions of the critical 
attractor. In practice, a partition of the phase space on Nt, ox cells of equal size is performed and a set of N c 
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identical copies of the system is followed whose initial conditions are uniformly spread over the phase-space. The 
ratio r = N c /Nb ox is a control parameter giving the degree of sampling of the phase-space. 

Within the non-extensive Tsallis statistics, there is a proper entropy S q evolving at a constant rate such that 

K q = lim [S q (t) - S q (0)]/t (5) 

N bom — too 

goes to a constant value as t — > oo. Notice that K q < for the process of convergence towards the critical attractor. 
Assuming that all cells of the partition are occupied with equal probability, the entropy S q (t) can be written as 

The last two equations imply that the number of occupied cells evolves in time as 

W(t) = [W{0) 1 - q + (1 - q^at} 1 '^ (7) 

with the exponent fi = — l/(q — 1) > governing the asymptotic power-law decay. 

In figure 1, we show our results for W(t)/Nb ox in the standard logistic map with inflexion z — 2 and from distinct 
partitions of the phase space with sampling ratio r = 0.1 . We observe that, after a short transient period when 
W(t) is nearly constant, a power-law contraction of the volume occupied by the ensemble sets up. W(t) saturates 
at a finite fraction corresponding to the phase space volume occupied by the critical attractor on a given finite 
partition. The saturation is postponed when a finer partition is used once the fraction occupied by the critical 
attractor vanishes in the limit Nb ox — y oo. 

In figure 2, we show W(t)/Nb ox for a given fine partition of the phase-space and distinct sampling ratios r. 
We notice that the crossover regime to the power-law scaling is quite short for large values of r so that a clear 
power-law scaling regime sets up even at early times. This feature is consistent with eq.(7) which states that the 
crossover time r scales as r ~ 1/W (Q) q ^ 1 . Further, the scaling regime exhibits log-periodic oscillations once the 
multifractal nature of the critical attractor is closely probed by such dense ensemble. A general form for W(t) 
reflecting the discrete scale invariance of the attractor can be written as 

where P is a function of period unity and A is the characteristic scaling factor between the periods of two consecutive 
oscillations. These log-periodic oscillations have been observed in a large number of systems exhibiting discrete 
scale invariance [Q. In general the amplitude of these oscillations ranges form 10~ 4 up to 1CP 1 . Keeping only 
the first term in a Fourier series of P(lni/lnA), one can write W(t) in the form 



W{t) = c t-» 



ci / hit 

2 — COS 2"7T- 



(9) 



co \ In A 

Log-periodic modulations correcting a pure power-law have been found in several systems, as, for example 



diffusion- limited-aggregation (T8| , crack growth [|T9| , earthquakes 1 20 1 and financial markets |2l]] . It has also been 
observed in thermodynamic systems with fractal-like energy spectrum j22 23 1. The factors controlling the log- 
periodic relative amplitude 2ci/co are not well known for most of the systems where it has been observed. In 
the present study, we can closely investigate the factors which may control these amplitudes by measuring it as 
a function of the map inflexion z for a fixed partition and sampling ratio (see figure 3). We found that these 
oscillations have amplitudes decaying exponentially with z as shown in figure 4. It is interesting to point out that 
the fractal dimension of the attractor is a monotonically decreasing function of z. Therefore, the above trend 
indicates a possible correlation between the amplitude of the log-periodic oscillations and the fractal dimension of 
the dynamical attractor. 

We also measured the critical exponent /i as a function of the map inflexion z. Our results are summarized 
in the Table. It is a decreasing function of z as can be seen in figure 5. The volume occupied by the ensemble 
depicts a fast contraction for z ~ 1 where the fractal dimension is small. On the other side, a very slow contraction 
is observed for large values of z, pointing towards a saturation or at most to a logarithmic decrease of W(t) in 
the limit of dense attractors. We would like to point out here that the exponent governing the expansion of the 
volume occupied by an ensemble of initial conditions concentrated around the inflexion point exhibits a reversed 
trend. Although scaling arguments have shown that this exponent can be written in terms of scaling exponents 
characterizing the extremal sets in the attractor, we could not devise a simple scaling relation between /i and the 
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multifractal singularity spectrum. However, we observed that, when plotted against the fractal dimension of the 
attractor as shown in figure 6, the dynamic exponent ^ is very well fitted by/itx(f — df) 2 , which indicates df 
as the relevant geometric exponent coupled to the dynamics of the uniform ensemble. We would like to mention 
here that the same dynamic exponents were obtained for the generalized periodic maps which belong to the same 
universality class of logistic-like maps |10|. 



III. THE CONVERGENCE TO THE CRITICAL ATTRACTOR OF THE CIRCLE MAP 

The results from the previous section indicate that a slow convergence to the critical attractor shall be expected 
for dense critical attractors. However, it is not clear in what fashion this convergence will take place when the 
dynamical attractor fills the phase space with a multifractal probability density as occurs for the one-dimensional 
critical circle map 

9 t+1 = d t + ft - ^- sin (27rf? t ) mod(l) (10) 

Z7T 

where < 9 t < 1 is a point on a circle. The circle map describes dynamical systems possessing a natural frequency 
uj\ which are driven by an external force of frequency loi (ft = uj\ / ll>2 is the bare winding number) and belongs to 
the same universality class of the forced Rayleigh-Benard convection For £1 = 0.606661... the circle map has 
a cubic inflexion (z — 3) in the vicinity of the point 9 = 0. Starting from a given point on the circle, it generates a 
quasi-periodic orbit which fills the phase-space and the dynamical attractor is a multifractal with fractal dimension 
dj 1 0- 

In figure 7 we show our results for the temporal evolution of the phase-space volume occupied by an ensemble of 
initial conditions uniformly spread over the circle. W(t) exhibits a rich pattern which resembles the one observed 
for the sensitivity function associated to the expansion of the phase-space from initial conditions concentrated 
around the inflexion point. However, W(t) does not present any power-law regime. Instead, the lower bounds 
display a slow logarithmic decrease with time, saturating at a finite volume fraction. The saturation is a feature 
related to the finite partition used in the numerical calculation. This minimum decreases logarithmically with 
the number of cells in the phase-space as shown in figure 8. We also observed the same behavior for generalized 
circle maps with an arbitrary inflexion z The critical attractors within this family have all d/ = 1 although 
they exhibit a ^-dependent multifractal singularity spectra. The ^-independent scenario for W(t) corroborates the 
conjecture that df is the relevant geometric exponent coupled to the dynamics of the uniform ensemble. 



IV. SUMMARY AND CONCLUSIONS 



In this work, we studied the temporal evolution in phase space of an ensemble of identical copies of one- 
dimensional nonlinear dissipative maps. We found that the phase-space volume occupied by an initially uniform 
ensemble displays a power-law decay with log-periodic oscillations whenever the dynamical attractor has a fractal 
dimension df < 1. The amplitude of the oscillations was found to depict a monotonic parametric dependence 
on df. For dense multifractal attractors, W(t) presents only a slow logarithmic contraction of its lower bounds 
followed by a rich pattern. 

The critical exponent characterizing the contraction of the uniform ensemble was found to have no direct relation 
to the one governing the expansion from a set of initial conditions concentrated around the inflexion point. In 
particular, no power-law was found for the contraction in the standard and generalized circle maps, in contrast 
to the ^-dependent power-law expansion. These results indicate that the relevant Tsallis entropy evolving at a 
constant rate (modulated by log-periodic oscillations) is characterized by an cntropic index q that depends on the 
initial ensemble. It would be valuable to investigate the possible existence of classes of ensembles with a common 
dynamics in phase-space and, therefore, characterized by the same entropic index q. The non-universality of q 
with respect to the initial ensemble is related to the multifractal character of the dynamical attractor. However, 
as for the ensemble concentrated at the vicinity of the inflexion point, the exponent governing the dynamics of 
the uniform ensemble is coupled to a geometric scaling exponent, in particular to the proper fractal dimension of 
the attractor. This result comes in favor of the concept that the degree of nonextensivity of the entropy measure 
evolving at a constant rate is related to the fractal nature of the dynamical attractor. 
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FIGURE AND TABLE CAPTIONS 



Figure 1 - The volume occupied by the ensemble W(t) as a function of time in the standard logistic map (z = 2) 
and with sampling ratio r = 0.1. From top to bottom N box = 2000, 8000, 32000, 128000. 

Figure 2 - The volume occupied by the ensemble W(t) as a function of time in the standard logistic map (z = 2) 
and for a partition containing Nb ox = 128000 cells. Notice the emergence of log-periodic oscillations for large 
sampling ratios. 

Figure 3 - The periodic function W (t) / '(cot~ A1 ) versus time within the scaling regime and for r = 10. Data 
from map inflexions z = 1.1, 1.25, 1.5,2.0 are shown. The amplitude of the oscillations decreases monotonically 
as z increases, but the characteristic scaling factor between the periods of two consecutive oscillations is roughly 
z-independent. 

Figure 4 - The amplitude of the log-periodic oscillations 2ci /co as a function of the map inflexion z for sampling 
ratio r = 10. The monotonic decrease of the oscillations indicates a close relation between these and the fractal 
dimension of the underlying dynamical attractor. 

Figure 5 - The dynamic exponent \x governing the contraction of the occupied phase space volume [W(t) oc 
as a function of the map inflexion z. 

Figure 6 - The parametric dependence of the dynamic exponent \x with the fractal dimension df of the critical 
attractor. It is very well fitted to the form fj, oc (1 — df) 2 , indicating that df seems to be the relevant geometric 
exponent coupled to the dynamics of the uniform ensemble. 

Figure 7 - The volume occupied by the ensemble W(t) as a function of time in the standard critical circle map. 
The lower bounds display a slow logarithmic decay with time saturating at a finite volume fraction due to the 
finite partition of the phase space. 

Figure 8 - The asymptotic lower bounds for the occupied volume in the phase space versus the number of cells 
Nb ox . The logarithmic decay agrees with the prediction that n(df — > 1) — ► 0. The same behavior was observed for 
the family of generalized circle maps and corroborates the conjecture that df is the relevant geometric exponent 
coupled to the dynamics of the uniform ensemble. 

Table - Numerical values, within the z-generalized family of logistic maps, of: i) the dynamic exponent /i 
governing the contraction towards the critical attractor of the uniform ensemble; ii) the entropic index q of the 
proper Tsallis entropy decreasing at a constant rate; Hi) the fractal dimension df of the critical attractor. These 
values also hold for the generalized periodic maps. The last line represents our results for the z-generalized circle 
maps. 
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Table 



z 


/i =-l/(l- g ) 


q 


df 


1.10 


1.62 ±0.02 


1.62 ±0.01 


0.32 ±0.02 


1.25 


1.23 ±0.01 


1.81 ±0.01 


0.40 ±0.01 


1.5 


0.95 ±0.01 


2.05 ±0.01 


0.47 ±0.01 


1.75 


0.80 ±0.01 


2.25 ±0.015 


0.51 ±0.01 


2.0 


0.71 ±0.01 


2.41 ±0.02 


0.54 ±0.01 


2.5 


0.59 ±0.01 


2.70 ±0.02 


0.58 ±0.01 


3.0 


0.515 ±0.005 


2.94 ±0.02 


0.60 ±0.01 


5.0 


0.395 ± 0.005 


3.53 ±0.03 


0.66 ±0.01 


z-circular 
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